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ABSTRACT 

The statistical distribution of the masses of planets about stars between the Sun 
and the center of the galaxy is constrained to within a factor of three by an intensive 
search for planets during microlensing events. Projected separations in terms of the lens 
Einstein ring radius yield a rough estimate of the distribution of planetary semimajor 
axes with planetary mass. The search consists of following ongoing stellar microlensing 
events involving sources in the center of the galaxy lensed by intervening stars with high 
time resolution, 1% photometry in two colors in an attempt to catch any short time 
scale planetary perturbations of the otherwise smooth light curve. It is assumed that 
3000 events are followed over an 8 year period, but with half of the lenses, those that are 
members of binary systems, devoid of planets. The remaining 1500 lenses have solar- 
system-like distributions of 4 or 5 planets. The expectations from the microlensing search 
are extremely assumption dependent with 56, 138, and 81 planets being detected for three 
sets of assumptions involving how the planetary masses and separations vary with lens 
mass. The events can be covered from 54% to 62% of the time on average by high time 
resolution photometry from a system of three or four dedicated two meter telescopes 
distributed in longitude, so 38% to 46% of the detectable small mass planets (very short 
perturbations of the light curve) will be missed. But perturbations comparable to a day 
in length means all of the detectable Jupiters and Saturns will in fact be detected as 
well as a large fraction of the Uranuses. The ground based observational technique is 
robust, and meaningful statistics on planetary masses and separations can be inferred 
from such an intensive search, although these statistics, like the inferred data set, will 
also be dependent on the assumptions about the nature of the set of planetary systems. 
Finding most of many giant and sub giant planets outside the Einstein ring radii of their 
respective stars may be a better indicator of the frequency of Earth mass planets than 
direct detection of a few of the latter. 
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1. Introduction 

The nearly circular and coplanar orbits of the plan- 
ets around the Sun point to the formation of these 
planets from a highly dissipative disk of gas and dust. 
The inevitability of the natural formation of such a 
disk during the gravitational collapse of a rotating 
molecular cloud to form a star and the observational 
confirmation that essentially all recently formed stars 
possess such a disk {e.g. Strom et al. 1995) have 
led most to believe that planetary systems are a com- 
mon result of the star formation process and therefore 
must be ubiquitous in the galaxy. The recent discov- 
ery of planetary sized bodies around several nearby 
stars (Mayor and Queloz, 1995; Marcy and Butler, 
1996; Butler and Marcy, 1996; Gatewood, 1996) and 
about the pulsar PSR 1257-1-12 (Wolszczan and Frail, 
1992) reinforces this confidence. However, except for 
the system around the M star Lalande 21185 (Gate- 
wood, 1996), the inferred planets are in systems very 
much unlike our own, since planets comparable to or 
exceeding Jupiter's mass orbit very close to their pri- 
mary stars. In fact, those companions in eccentric 
orbits about the ordinary stars are inferred to not be 
planets at all but were formed like binaries by frag- 
mentation of a collapsing molecular cloud (Boss, 1996; 
Mazeh et al. 1996). This leaves massive planets com- 
parable to or exceeding Jupiter's mass in nearly cir- 
cular orbits about 51 Peg (msini ~ 0.45mj, 4.2 day 
period), 55 Cnc (msini ~ 0.78m.7, 14.7 day period) 
and 47 UMa (msini ~ 2.4mj, 1090 day period) con- 
trary to our expectations of the stellar distances at 
which giant planets would be formed. (However, there 
is a mechanism by which such planets could migrate 
toward their stars after formation provided a nebula 
with significant mass persists sufficiently long (Lin, 
et al. 1996; Ward and Hourigan, 1989; Ward, 1996)). 
These discoveries imply that other planetary systems 
are likely to differ greatly from our own, with totally 
different distributions of planetary mass and charac- 
ter with distance from the star. In particular the 
number of terrestrial type planets in any system may 
be considerably different from expectations based on 
our own solar system. 

The theory of the formation and evolution of plan- 
etary systems depends on only one example, whereas 
the recent discoveries of other planets indicate that a 
variety of scenarios leading to systems strikingly dif- 
ferent from our own may be common. To understand 
the frequency of the particular chain of events that led 



to a planetary system like ours, where life was able to 
develop on one terrestrial planet, we need statistics 
not only on the occurrence of other planetary sys- 
tems, but also on the planetary mass distributions as 
a function of the planet separation from the parent 
star. Giant planets close to their primaries may pre- 
clude the existence of terrestrial type planets in those 
particular systems. The recent discoveries of planets 
around other stars depended on detecting motion of 
the star relative to a system center of mass with the 
period of the planetary orbit. High precision Doppler 
spectroscopy and, in the case of Lalande 21185, high 
precision astromctry were used to determine the stel- 
lar motions. These techniques can provide informa- 
tion on planetary mass distributions and separations, 
where the Doppler technique is most sensitive to plan- 
ets that are close to their stars and astrometry most 
sensitive to planets that are far away. Unfortunately, 
these observations arc very demanding so relatively 
few stars can be monitored, and the observations must 
be continued for at least a planetary period for secure 
detection. Several decades of dedicated effort would 
be required to compile adequate statistics, and nei- 
ther method is sensitive to terrestrial- mass planets at 
any separation. 

The recently announced intention of the National 
Aeronautics and Space Administration to search for 
planets about other stars (Beichman, 1996) has as 
its goal the detection and characterization of terres- 
trial type planets orbiting stars within about 10 pc 
from the Sun. The instrumentation to be developed 
in this search depends on the distribution of planetary 
masses and in particular on the frequency of occur- 
rence of Earth-mass planets. Planning this program 
requires the planetary statistics on a relatively short 
time scale, and it requires information on the fre- 
quency of occurrence of terrestrial-mass planets that 
is not obtainable by the radial velocity and astrom- 
etry techniques. Mao and Paczyhski (1991) pointed 
out that a planet could be detected as a companion of 
a star by perturbing the light amplification curve of a 
more distant star (source) that is being gravitation- 
ally lensed (microlensing) by the nearer star (lens). 
Later Gould and Loeb (1992) derived a probability of 
nearly 20% for detecting a Jupiter about a solar mass 
star given that a microlensing event was taking place 
by assuming that a 5% perturbation of the light curve 
was detectable. Microlensing is a rare event, so mil- 
lions of stars must be monitored in order to catch the 
few whose apparent brightnesses change in a system- 



2 



atic, achromatic way due to microlensing by a nearer 
star. The MACHO (MAssive Compact Halo Objects) 
coUaboration has been remarkably successful in de- 
tecting approximately 100 microlensing events toward 
the galactic bulge over a two year period. The pos- 
sibility of detecting at least 350 events per year with 
modest upgrades in technology together with a larger 
dedicated telescope at a good site (C. Stubbs, 1996, 
private communication) promises a means to gather 
planetary statistics at a reasonably rapid rate. 

Although planets that orbit very close to their stars 
will not be detectable during a microlensing event, 
(The star and planet would act as a single lens.), 
the likely situation that most planetary systems will 
be dispersed more like our own motivates consider- 
ation of a dedicated microlensing search for planets. 
How many planets are likely to be detected in such a 
search, what is the nature of the data set, and what 
planetary statistics can be obtained? Here we con- 
struct examples of answers to these questions based 
on the Gould and Loeb (1992) probability of detecting 
a planet given that the planet's central star is acting 
as a lens during a microlensing event. 

The source stars are treated as point sources in 
the derivation of this probability, which means that 
the application to small mass planets is not really 
valid. However, Bennett and Rhie (1996) have calcu- 
lated detection probabilities given that a microlensing 
event is taking place for planet/ (stellar lens) mass ra- 
tios m/M = 10~^ and lO"** at distances of 4 and 6.4 
kpc and for sources at the galactic center of radii 3i?0 
and 13i?0 as a function of the planet-star separation. 
These detection probabilities averaged over the mass 
function of the lenses for the 3i?Q source radius and 
m/M = 10~^ for a set of specific planetary orbital 
semimajor axes are remarkably close to the proba- 
bilities for the same fixed mass ratio from the point 
source calculation (Table ||) . This indicates that our 
probabilities of detection for the Earth-mass planets 
for point sources may be representative to within a 
factor of 2. Even without this near agreement, our 
purpose here is to demonstrate the assumption de- 
pendence of the expectations from any microlensing 
search for planets, to show the effect of the stellar 
mass function where most of the lenses will be spec- 
tral class M stars, and to show how well a minimum 
data set can be interpreted. The variations in the 
number of planetary detections due to a variation in 
the assumptions is at least comparable to the errors 
in the estimates due to considering only the conse- 



quences for point sources. 

We begin in Section 2 with a brief description of 
the simplest microlensing theory, and a description of 
the Gould and Loeb (1992) derivation of the proba- 
bility of detecting a planet during a Icnsing event suf- 
ficient to allow one to understand the scaling of the 
published probability to arbitrary planet/star mass 
ratios and arbitrary mean lens distances. The distri- 
bution of events over a lens mass function like that 
of the stars in the local neighborhood is developed in 
Section 3, and it is used to average the scaled prob- 
ability of detection of a planet of given mass, orbit 
semimajor axis, and average distance of the lens from 
the observer. 

A sample observing program like that advocated 
in the Roadmap for the Exploration of Neighboring 
Planetary Systems (Tytler, 1995) is described in Sec- 
tion 4 after showing how the planet/star mass ra- 
tio and the planet-star separation normalized by the 
Einstein ring radius of the lens star is deduced from 
the microlensing light curve perturbed by the planet. 
Two color, high time resolution photometry of all on- 
going microlensing events with 1% accuracy is ad- 
vocated in this program. Two colors are necessary 
for redundancy, for distinguishing intrinsically vary- 
ing stars, and to reduce uncertainties in some cases 
where the source star is resolved by the planet. The 
probability of being able to observe the high density 
star region in the galactic center (GC) from each of 
four specific observing sites at any specific time is used 
to develop the overall probability of observing the GC 
from at least one of the observatories as a function of 
the time of year. This is averaged over the time the 
GC is observable to show that this particular distri- 
bution of telescopes for a microlensing search could 
cover the ongoing events with high time resolution 
photometry about 54% to 62% of the time. 

A reasonable data base of 3000 events assumed col- 
lected over an 8 year period is used in Section 5 to de- 
termine the number of planets that would be detected 
for three sets of assumptions about systems similar to 
the solar system. For one of these systems a sample 
planetary data set consisting of the planet/star mass 
ratios and planet-star separations is constructed from 
the appropriate probability distributions. An inter- 
pretation of this data set follows in Section 6 where 
it is assumed that only information from the detailed 
light curve is available with spectral classification of 
the brighter lenses used to further constrain these lens 
masses. Although only the planet/star mass ratios 
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follow from the light curves, the limited range of lens 
masses means a reasonable estimate of the distribu- 
tion of planetary masses is obtained. Observational 
selection effects make the orbital semimajor axes of 
the planets less defined except for the general trend 
relative to the Einstein ring radius that reflects the 
assumptions about the nature of the systems. 

In Section 7 we show how the finite sizes of the stel- 
lar sources affects the detection of Earth sized plan- 
ets, where one gains because of the higher probabil- 
ity of the small planet affecting the light curve and 
longer duration of the event, but loses because of the 
much lower amplitude of the event. We also show 
examples of averaging the Bennett and Rhie (1996) 
detection probabilities (for finite sized sources) over 
the lens mass function, compare these probabilities 
to those obtained for point sources and point out the 
necessary steps for obtaining an overall probability of 
detecting small mass planets while accounting for the 
finite source sizes of the galactic bulge stars. 

2. Microlensing and planet detection 

The term microlensing comes from extra-galactic 
astronomy where the images of a quasar lensed by 
a single star in an intervening galaxy are separated 
by the order of microarcseconds (Chang and Refsdal, 
1979). The images of a star in the galactic bulge that 
is lensed by an intervening star are separated by mil- 
liarcseconds, but the term is still applied. Fig. |l| 
shows the geometry of the gravitational lensing of a 
point source and the definition of several parameters. 
O, L and S refer to observer, lens and source respec- 
tively. The angle a is the general relativistic (small 
angle) deflection of a ray from the source that reaches 
the observer. The point source lens equation deduced 
from the figure is — To r — i?| 0, where 

R. ^ (1) 

is the Einstein ring radius. The Einstein ring is the 
symmetric image of the source when it is directly be- 
hind the lens (vg — 0). Here G is the gravitational 
constant, M is the mass of the lens, c is the velocity 
of light and z = Dq^/Dqs. The circle centered on L 
in Fig. is the Einstein ring on the lens plane. So- 
lution of the lens equation ri.2 = (ro ± \/rl + 4i?|)/2 
yields the positions of the two images Ii and I2 out- 
side and inside the Einstein ring respectively. The 



amplification {e.g., Gould and Loeb, 1992), 

is the factor by which the flux density from the source 
star is increased. It is found for each image as the 
transformation of areas from source to image coordi- 
nates with the results for both images being added 
together. Here u = Tq/Re. The light curve of the 
event reaches a maximum when u has its minimum 
value of Ufc. The time scale of the event is the time 
for the projected position of the source to move an 
Einstein ring diameter, r = 2Re/v oc ^/M where v 
is the relative speed of the source projected onto the 
lens plane. The singularity in A at u = is called the 
caustic point. This singularity is removed when the 
source is not a point source. Fig. || shows the light 
curves for a variety of impact parameters and illus- 
trates the independence of the perceived time scale 
on the latter. 

If the lens is a relatively close binary system, the 
lens equation relating coordinates of the images in 
the lens plane to the coordinates of the source in the 
source plane is now a vector equation (Schneider and 
Weiss, 1986). The quadratic lens equation for a single 
lens becomes two 5th order equations or a single fifth 
order equation in the complex plane (Witt, 1990). 
There are now either three or five images and the 
amplification is calculated for each image as before as 
the transformation of areas with the results summed 
for the total. The caustic point for the single lens 
is now one or two closed caustic curves in the source 
plane that transform into critical curves in the lens 
plane, where A — > 00 when the point source is on a 
caustic. The singularity again vanishes for finite sized 
sources. There are three images when the source is 
outside a closed caustic curve and five when it is in- 
side. A sharp peak in the light curve occurs when the 
source crosses a caustic with increased amplification 
while the source is inside. The simple, bell shaped 
light curve for a single lens can become quite com- 
plicated for a binary lens but also quite informative 
{e.g. Alard, et al. 1995). 

When the second member of the binary system is 
a planet, the small mass of the planet compared to 
the star leads to a significant deviation of the light 
curve from that of the star without the planet only if 
the one of the two unperturbed images of the source 
comes close to the planet's projected position in the 
lens plane (Gould and Loeb, 1992). Fig. || shows an 
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Fig. 1. — Geometry of a microlcnsing event and defi- 
nition of parameters. The circle in the lower diagram 
is the Einstein ring of the lensing star and Ii , and I2 
are the images of the source S. O, L and S refer to 
Observer, Lens and Source respectively. 
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Fig. 2. — Microlcnsing light curves for several impact 
parameters of the relative motion of source and lens 
stars. The independence of the time scale of the event 
on the impact parameter is illustrated. 



Fig. 3. — Example of a microlcnsing light curve 
perturbed by a planet. The planet /star mass ratio 
m/M — 0.001 and the planet is located 1.3i?£; from 
the lens. (After Gould and Loeb, 1992) 

example of such a perturbed light curve, where the 
planetary signature is larger than the simple sum of 
the gravitational lensing of the star plus the planet 
because the planet is focusing light already focused 
by the central star. For an event involving a lens 
with a planetary companion, the planet can affect 
the light curve of a point source by at most a very 
short fraction of the total time of the event as shown 
in Fig. 1^, and for most events, the planet will not 
reveal its presence. However, for a point source near 
the galactic center being lensed by a solar mass star 
at 4 kpc from the Sun, the probability that a Jupiter 
mass planet 5.2 AU from the star will cause at least a 
5% deviation from the unperturbed light curve some 
time during the event is about 0.17 (Gould and Loeb, 
1992). (An event is said to occur when the lens and 
source come within an Einstein ring radius of the lens 
of each other.) This high probability results from 
the semimajor axis being near the Einstein ring ra- 
dius of the star where the probability is maximized, 
but there is a considerable range of semimajor axes 
{2.4 AU < a < 7.3 AU for this case) where the prob- 
ability exceeds 0.1 (Gould and Loeb, 1992). David 
Tytler (private communication, 1995) has coined the 
term "lensing zone" for this range of semimajor axes 
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Fig. 4. — Analytic representation of the probability 
of detecting a planet with m/M — 0.001 given that 
an event is taking place, where a point source is as- 
sumed and where a 5% perturbation of the light curve 
anytime during the event is assumed detectable. The 
curve has been linearly extended beyond the Gould 
and Loeb data. The second form of is to account 
for the scaling of the position of the Icnsing zone for 
various average fractional distances of the lens Zg and 
lens masses M. {M is normalized by Mq.) 

where the probability of detection is high. 

An analytic representation of the Gould and Loeb 
(1992, Fig. 8) detection probability for a planet/star 
mass ratio m/M = 0.001 and for a 5% detection 
threshold (Appendix A) is shown in Fig. ^, where 
= ac/yJGM Dos is the scaled semimajor axis, 
and where the curve has been linearly extended for 
(3 < a;* < 5). We shall use this probabihty of detec- 
tion in generating a sample data set for a microlens- 
ing planet detection program, but several qualifica- 
tions and extensions are necessary. First the curve 
was derived under the assumption that the sources 
are point sources. The probabilities of detection are 
thus representative only as long as the angular size of 
the Einstein ring of the planet is large compared with 
the size of the source. This follows because a point 
source must be inside the Einstein ring of a lens for 
significant amplification and therefore only a fraction 
of a finite source can be amplified by a small lens 



with a consequently smaller perturbation of the light 
curve. At the same time, there is a larger probabil- 
ity that some part of the now finite sized image of 
the source will encounter the planet in the lens plane. 
Recent calculations by Bennett and Rhie (1996) have 
demonstrated that this enhanced probability of planet 
encounter may more than compensate for the lower 
amplitude of the perturbation in the sense that their 
maximum probability of detection is about twice the 
maximum probability of detection for a point source 
with a planet/star mass ratio of 10~^. Averaging the 
Bennett and Rhie detection probabilities over the lens 
mass function in Section 7 yields probabilities of de- 
tection close to those obtained with a point source, 
so small mass planet detections determined with the 
point source assumption should be within a factor of 
2 of those accounting for the finite source size. 

Next, the Gould and Loeb probability was derived 
by averaging over a distribution of lenses along a par- 
ticular line of sight that was displaced about 4° from 
the direction to the galactic center yielding a uniform 
distribution of lenses. Other lines of sight for the as- 
sumed spatial distribution of stars in the galaxy would 
yield different probabilities, but not significantly so. 
The fact that the number of detected events toward 
the galactic center was about three times that ex- 
pected from the distribution of stars assumed in the 
Gould and Loeb analysis (Alcock, et al. 1995) will 
affect the probabilities of detection by moving the po- 
sition of the average lens closer to the center of the 
galaxy thereby reducing the size of the Einstein ring 
radii (Eq. (|l|)). This follows since the excess lenses 
are thought to be located in a bar-like distortion of 
the galactic bulge distribution of stars oriented to- 
wards the Sun (Zhao et al. 1995). An approximate 
correction for this effect is made below. 

The correction in the probability for moving the 
average lens closer to the source as well as scaling the 
probability for different values of m/M will be better 
understood if we outline how the probability of de- 
tection is determined. If m/M 0, the binary lens 
solutions go to the single lens case with the two un- 
perturbed images described above. For small m/M 
only one of the images will be significantly affected 
when the planet approaches a separation from the 
unperturbed image position comparable to planet's 
own Einstein ring radius, and the new images will be 
close to the unperturbed position of the image. The 
amplification is calculated as before as the sum of the 
amplifications of each image. As the planet must be 



6 



close to the unperturbed image position for significant 
perturbation of the Hght curve, contours of planetary 
positions for a given percentage perturbation of the 
light curve will be grouped close to the unperturbed 
image position with the percentage perturbation de- 
creasing for contours further away from that position. 

In Fig. ^ we can imagine a contour of the plan- 
etary position in the lens plane corresponding to a 
5% perturbation surrounding each of the two nonper- 
turbed image positions. Examples of these contours 
are shown in Fig. 2 of Gould and Loeb (1992). If the 
planet is inside one of these contours drawn in the 
lens plane, the perturbation of the light curve will 
exceed 5%. As the source and lens pass, the unper- 
turbed images describe arcs in the lens plane (see Fig. 
0), and the loci of the contours as the images change 
their positions define two areas in the lens plane. If 
the projected position of the planet is inside either of 
these areas, there will be a perturbation of the light 
curve exceeding 5% sometime during the event and 
the planet will be detectable. The probability of the 
planet detection sometime during the event is just 
the probability that the planet is within either of the 
two areas. Finally the probability is determined for 
and averaged over all impact parameters less than the 
Einstein ring radius. 

The first step in determining this probability is to 
assume that a planet with semimajor axis a has a uni- 
form probability distribution over a sphere of radius 
a to account for random orbit inclination and phase. 
Then the probability that the planet has a projected 
separation between r and r + dr is 



probability that an event will occur to obtain 



p(r)dr 



rdr 



(3) 



and the probability that the planet is between r and a 
is F{r) = y/l — ja?' ■ If we normalize the distances 
by the Einstein ring radius Re from Eq. (|^), 



l-4z(l-z)^ 



1/2 



(4) 



represents the probability that the planet will be lo- 
cated between = r/RE and a/ Re from a lens 
located at fractional distance z to the source with 



= ac/ ^GMDqs as before. Gould and Loeb (1992) 
average Eq. (0) over the line of sight weighted by the 



F{xp) = 



/q^^ ^° ■■■dz + RE{z)v{z)F{xp,z)dz 

RE{z)v{z)dz 

(5) 

where i^(z) is the spatial density distribution and 
where Zo = if a;p < and Zg — 0.5^1 — xl/x^ 

if Xp > x^. Although this is the form used by Gould 
and Loeb in their calculations, it differs from their Eq. 
(3.8) in that the integral is separated into two seg- 
ments. The split limits are necessary since the inte- 
grand becomes undefined when Xp > x^. At z = 1/2, 
Xp > x^ implies r > a which is clearly impossible. 
As z moves away from 1/2, Xp — t/Re gets larger as 
Re gets smaller. So to have a value of r still inside 
the maximum planetary separation from the star for 
larger Xp, the integration can extend only over those 
values of z where large Xp still corresponds to r < a. 

The particular form of ^{z) chosen by Gould and 
Loeb places the average lens near 4 kpc from the Sun. 
If we let v{z) = C5{z ~ z^) where C is some constant, 
all of the lenses are at z = Zs and F(xp) — F(xp, Zs) = 
yj\ — 4zs(l — Zs)(xp/x,)^. The probability distribu- 
tion 



dF(xp) 



^XpZs{\ - Zs) 



^Zs{l - Zs){Xp/x^y 



is strongly peaked at a removable singularity at a;, = 
Xp \J 4zs(l — Zs). This is where the probability of find- 
ing the planet per unit area in the lens plane would 
be maximal. Now the contours surrounding the un- 
perturbed image positions within which the perturba- 
tion to the light curve exceeds 5% are maximal in size 
when the image comes close the the Einstein ring and 
the locus of contour areas during the event is likewise 
maximized. As the planet must be near the image 
to be inside the contours, Xp ^ \ (i.e. planet near 
the Einstein ring). The combination of the largest 
areas in the lens plane for which the perturbation of 
the light curve would exceed 5% sometime during the 
event, if the planet were inside the area, with the 
maximum probability per unit area corresponds to 
the highest probabilities for finding the planet. So for 
Xp near 1, ~ \J 4zs(l — z^) corresponds to the peak 
in the probability curve. For z^ = 1/2 the peak should 
be near x» = 1, and we see from Fig. ^ that placing 
all the lenses at z = 1/2 gives a peak in the probabil- 
ity of detection near the peak that was obtained for 
the average lens being at z = 1/2. So for the average 
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lens being at a value of z other than 1/2, we shift 
the lensing zone by \/ 4z(l — z). This is equivalent 
to writing a;* = 0.25a/ y^4Mzs(l — Zg) in reading the 
probabilities from Fig. |[ where a is in AU, M is in 
units of Mq , Zj here is approximately the value of z 
for the average lens, and other constants have been 
evaluated. 

Next we must scale the probability of detection to 
planet/star mass ratios m/M other than the value of 
0.001 appropriate to Fig. ^. We note that the size of 
the 5% contour around the unperturbed image posi- 
tion in the lens plane scales as m/M (Gould and 
Loeb, 1992). As the arc described by unperturbed 
images is long compared to any dimension across the 
contours and this arc is independent of m/Af, only 
the dimension of the area of the locus of contours per- 
pendicular to the image arc is affected by a change in 
m/M, so the area of the locus of 5% contours for the 
complete event scales as yJrnjM. Since the proba- 
bility per unit area of the planet falling into the area 
of the locus of contours does not vary much over the 
area, the probability is changed approximately as the 
area is changed. Hence, for a given value of x*, we 
need but multiply the probability by ^/"i/O.OOlM to 
obtain the probability of detection for an arbitrary 
planet/star mass ratio. We will thus represent the 
Gould and Loeb probabihty curve in Fig. || so scaled 
and shifted as P{a,m, M, Zs). This is the probabil- 
ity, during a lensing event, of detecting a planet of 
mass m, with semimajor axis a, orbiting a lens star 
of mass M located at a mean fractional distance z^. 
Below we will average this probability over the mass 
function for particular sets of values of a,m, z^, but 
first the distribution of events over the mass function 
must be developed. 

3. Distribution of the events over the lens 
mass function 

We shall assume that the mass function for local 
stars is representative of the stars between us and 
the center of the galaxy. This assumption will prob- 
ably be most in error for masses larger than about 
1.5Mq, since more massive lenses within the galac- 
tic bulge will have evolved off the main sequence and 
thereby will have lost mass on their way to becom- 
ing white dwarfs or neutron stars. The mass function 
used here is an analytic approximation derived from 
the multiplicity corrected present day mass function 
for main sequence stars given in Table 1 and Fig. 5 of 



Basu and Rana (1992). We consider only the range 
of masses 0.08 < AI < 2.0 because of the very small 
fraction of stars of higher mass and because of the 
completed evolution of the older possible lens stars in 
the near side of the galactic bulge. M is expressed 
in units of the solar mass Mq . This mass function of 
main sequence stars neglects the contribution of white 
dwarfs and neutron stars to the lens population, but 
these latter stars are expected to be a small fraction 
of those stars remaining on the lower main sequence. 
The mass function is represented by (based on Basu 
and Rena, 1992) 

(t){M) = 48.39Af"i 0.08 < M < 0.5327, 
= 13.73M"3 0.5327 < Af < 1.205, 
= 28.04A/"'5-^^ 1.205 < M < 2.0, (6) 

where (j){M)dM is the number of stars/pc^ in the solar 
neighborhood in the mass range dM about M inte- 
grated perpendicular to the galactic plane. 

The fraction of stars in the range dM about M 
is found by dividing (j>{M)dM by the integral of 
4>{M)dM over 0.08 < M < 2.0, but this is not the 
fraction of microlensing events to be expected for 
lenses in mass increment dM. This follows from our 
definition of an "event" as the lens and source com- 
ing within an angular radius of each other equal to 
the Einstein ring angular radius Re/Dql of the lens. 
The dependence of on M (Eq. @)) means that a 
more massive lens has a greater probability of having 
a source pass within its Einstein ring radius than a 
less massive lens. 

To evaluate the distribution of events over the 
mass function, we assume first for simplicity that 
Dos — 8kpc is the same for all the source stars in 
the galactic bulge, so that z = Dol/Dqs is the frac- 
tional distance of the lens to the center of the galaxy. 
If (z) is the magnitude of the component of velocity 
of a single lens perpendicular to the line of sight, and 
N is the number of sources per unit solid angle toward 
the galactic bulge, Nv^{z)2RE{M,z)/{z'^Dlg) is the 
number of encounters (events) per unit time experi- 
enced by the single lens of mass M at fractional dis- 
tance z. Next, let v{r) be the number density of stars 
at distance r, such that u{r)dr is the areal density of 
stars for the range between r and r+dr. Let f{M)dM 
be the fraction of stellar masses between masses M 
and M + dM, which is constructed from (t>{M) (Eq. 
(^) as described above. Setting r = DqsZ, we have 
f{AI)dMi'{z)DosZ^dz is the number of stars between 
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M and M + dM per unit solid angle that are in the 
slab between z and z + dz. Multiplying this by the 
events per unit time due to a single lens of mass M 
at z and inserting the explicit M and z dependence 
into Re yields 



^^DosNf{M)dM 



v{z){v^{z)) V z(l - z)dz 
(7) 

as the number of events per unit solid angle per unit 
time due to lenses between M and M + dM over the 
entire line of sight. Here {vj) is an average transverse 
speed. Then from Eq.(|^), 



T{M)dM 



Mf{M)dM 

1.08 



Io:os^fiM)dM 



(8) 



represents the fraction of events due to lenses with 
masses between M and M + dM, where common fac- 
tors from Eq. (^) have been cancelled in numerator 
and denominator. The probability of detection of a 
planet of mass m, semimajor axis a, located at aver- 
age fractional distance Zg weighted by the distribution 
of events over the mass function along the line of sight 
is 



P{m,a,Zs) 



T{M)P{m,M,a,Zs)dM. (9) 



0.08 



The distribution of this probability over the m, a 
plane is shown in Fig. ^ for the case where the plane- 
tary mass m and semimajor axis a are invariant over 
the lens mass function. It is used in section 5 to 
deduce sample results from a dedicated microlensing 
search for planets, but first in the following section 
we define a representative observing program that ac- 
counts for observing constraints. 

4. Sample Observing Program 

Recall that the planet-image separation must be 
on the order of the Einstein ring radius of the planet 
or less to cause a significant perturbation of the sin- 
gle lens light curve. The time scale for the whole 
event is the time required by the source to traverse 
an Einstein ring diameter of the star (see Fig. The 
time scale of the planetary perturbation is the time 
required for the image to traverse the Einstein ring di- 
ameter of the planet, which is comparable to the time 
required for the source to traverse the Einstein ring 
diameter of the planet. The relative transverse speed 



of the source is the same for both cases, so the ra- 
tio of the time scale of the planet perturbation to the 
time scale of the event is just ^Jm/M (see Eq. If 
an event perturbed by a planet with mass ratio 0.001 
is observed to last 45 days, the planet perturbation 
would last only about 1.4 days. For a point source, a 
perturbation by an Earth mass planet about a solar 
type star (m/M = 3 x 10~^) would last only about 
1.8 hours for the same 45 day main event. However, 
the effect of the finite size of the source on the de- 
tection of small mass planets is profound, where a 
perturbation of the light curve is both reduced in am- 
plitude and extended in duration beyond the 2 hours 
expected for a point source. For this reason, the ob- 
serving strategy must be fundamentally different for 
detecting the small mass planets. We will discuss this 
in more detail in Section 7. If a 5% perturbation of a 
light curve is considered detectable, all events would 
have to be photometrically monitored with about 1% 
accuracy at time intervals no longer than about one 
hour to detect a perturbation by Earth-mass planets 
and if a perturbation is detected, at intervals of a few 
minutes thereafter on that particular light curve for 
the duration of the perturbation. In reality, the time 
scale for the planetary perturbation for a given mass 
ratio will vary somewhat depending on the angle that 
the source crosses the given magnification contours 
in the source plane. But that angle is known from 
the position of the planetary perturbation on the un- 
perturbed light curve. This defines the intersection 
of the source trajectory with the line connecting the 
lens and planet. From this angle and an accurate, 
high time resolution light curve obtained during the 
perturbation, one can model the light curve and refine 
the determination of the mass ratio. 

The ratio of the planet-star separation and the Ein- 
stein ring radius rp/R^ also follows from the posi- 
tion of the planetary perturbation on the unperturbed 
light curve. If the perturbation is removed, the unper- 
turbed amplification yields u — r/R^ from Eq. (H). 
For each value of u, the positions of the two unper- 
turbed images are known and the planet must be near 
one or the other image to have perturbed the light 
curve. The selection of the correct image between the 
two can follow from the modeling of the light curve. 
If the caustic in the source plane is intercepted by 
the point source, the planetary position is known to 
within the Einstein ring radius of the planet. 

David Tytler (1995) has suggested a telescope 
distribution and observing program to accomplish 
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Fig. 5. — Probability of a usable night as a function 
of the time of year for four observatories that are suit- 
able locations for dedicated telescopes in a microlens- 
ing search for planets. 

the necessary observations. One of four two meter 
telescopes at established sites distributed as widely 
as possible in longitude in the southern hemisphere 
would monitor several tens of millions of stars in the 
galactic center (GC) in two colors at least once a night 
whenever the bulge is observable. This "survey tele- 
scope" would discover as many microlensing events 
as possible using a technique that is currently being 
used successfully by the MACHO collaboration (Al- 
cock et al. 1995, 1996). The three "foUowup tele- 
scopes" would photometrically monitor each of the 
ongoing microlensing events in two colors at the high 
accuracy and high time resolution described above. 
We have arbitrarily added a fourth foUowup telescope 
in Hawaii. 

The locations of three developed observatory sites 
in the southern hemisphere separated maximally in 
longitude plus Hawaii in the northern hemisphere are 
shown in Table 0. Weather, equipment problems 
and maintenance prevent 100% coverage of the events 
with high precision, high time resolution photometry. 
The fraction of usable nights (spectroscopic) at each 
observatory site is shown in Fig. H where the curves 
are constructed from monthly averages for LaSilla 
(http://lwhl.ls.eso.org) and South Africa (J. Cald- 



well, Private communication, 1996), quarterly aver- 
ages from Hawaii (J. Glaspey, Private communica- 
tion, 1996), and yearly averages from Siding Springs 
where summer and winter are comparable (J. Mould, 
Private communication, 1996). The GC can only be 
observed from a given site at some time t when it 
has a zenith angle less than 70° while the solar zenith 
angle is greater than 105°. If these criteria are satis- 
fied, the probability that the GC is observable at the 
particular time is taken from Fig. |^. Otherwise the 
probability of observing the galactic center from that 
site at that time is zero. 

The probability that the GC is observable from at 
least one of the observatories at any instant is then 
just the union {e.g. Feller, 1957, p. 88) of the prob- 
abilities for each site at that instant. This union of 
probabilities differs from that for one or another of the 
sites only when the GC is simultaneously observable 
from more than one site. The sum X^iLi ^(^j)^^i 
gives the fraction of a day that the GC can be ob- 
served, where P{ti) is the probability of observation 
at time ti and At^ = 1/n is an increment of a day with 
n being large but otherwise arbitrary. The probabil- 
ity of being able to observe the galactic center on a 
given day so calculated and thus to make photometric 
measurements of lensed stars is shown in Fig. ^ as a 
function of time. The maximum observability occurs 
in the southern winter and drops to zero when the 
Sun is near the position of the galactic center during 
southern summer. There would be a period of time 
spanning day 90 in Fig. ^ where the probability of ob- 
serving the GC would be unity if the weather were al- 
ways perfect and equipment always worked, since the 
GC would always be observable from at least one of 
the sites throughout the day. The addition of Hawaii 
to the three southern sites gives a significant increase 
in the probability of observation because of its low 
northern latitude and because of its ideal location in 
longitude between Siding Springs and LaSilla. 

As the Sun approaches the position of the galactic 
bulge in the sky, the duration of dark sky while the 
GC is up shrinks to the point where all of the ongo- 
ing events cannot be covered with high precision pho- 
tometry in the time available. In fact there is a point 
beyond which the survey telescope cannot cover all 
of the fields in its program, and observations would 
probably cease for a time longer than the 71 days 
when the GC is completely blocked by the Sun in 
Fig. ^. The MACHO group requires about an hour 
to cover 24 fields toward the galactic bulge with an 
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Observatory Latitude (deg.) E. Longitude (deg.) 
LaSilla ^29^25 289.27 
Siding Springs -3L27 149.06 
South Africa -32.38 20.80 
Hawaii +19.83 204.52 



Table 1: Location of selected observatory sites. 




100 200 300 

Time (Days from vernal equinox) 



Fig. 6. — Probability of observing the galactic center 
(GC) on a given day. A probability of 1 would indi- 
cate that the GC was observable at all times during 
the day from one or more observatories. The GC is 
considered observable if the zenith angle of the GC 
is less than 70°, the zenith angle of the Sun is more 
than 105° and the night is usable. 

additional hour of overhead (Alcock, et al. 1966; C. 
Alcock, private communication, 1996). If we assume 
that two hours is the minumum amount of dark time 
for observing the GC, the 71 day exclusion region in 
Fig. Ijis expanded to about 120 days, leaving 245 days 
as the maximum span of time per year for observing 
the GC in a microlensing survey and for foUowup pho- 
tometry. If we center the 245 day span over the peak 
probability in Fig. ^ and average the probability over 
this span of time, we obtain the average fraction of the 
time when the galactic center can be monitored with 
high time resolution photometry during observed mi- 



crolensing events. This fraction of 0.54 for just the 
three southern observatories is increased to 0.62 with 
the addition of a fourth telescope in Hawaii. 

The approximately 58% average coverage for the 
assumed configuration of telescopes means that a lit- 
tle less than half of the planetary perturbations of mi- 
crolensing light curves of very short duration will be 
missed. However, if the duration of the perturbation 
to the light curve exceeds a day, as it would in most 
cases for a Jupiter mass planet, the curves in Fig. || 
represent the fraction of the perturbation that could 
be monitored on the average where the probability of 
detecting the perturbation is essentially unity. This 
latter assertion follows since the probability of all of 
the sites having unusable nights on a given day is 

Pallcloudy = -P^)^3X IQ-HoS X IQ-*, 

i 

where Pi is the probability of a usable night at the ith 
observatory taken from Fig. |^, and where the range 
of values results from the variation in the weather 
conditions during the year. The fraction of the longer 
perturbations of the light curve that must be followed 
for meaningful interpretation is not large if the groups 
of high precision points are well distributed as they 
normally would be as the event is passed from one 
observatory to the next. 

A few more of the small planets will be missed be- 
cause of abandoning complete coverage to take mea- 
surements of ongoing perturbations every few min- 
utes, but these will be neglected in the following. The 
probability of observing the GC in Fig. || would ap- 
proach unity for a time spanning the southern winter 
if an arbitrary number of followup telescopes could be 
distributed more or less uniformly in longitude near 
-30° latitude {e.g., see Gould and Loeb, 1992). The 
impracticality of developing many new sites and nec- 
essary ocean gaps preclude this ideal. 

Alcock et. al. (1996) found 43 likely microlens- 
ing events in the first year galactic bulge data where 
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12.6 million stars were monitored for 190 days. For 
the same detection efficiency, following a little more 
than 100 million stars would allow the detection of 
350 events per year. However, from the Alcock et al. 
data, the optical depth for microlensing events for the 
distance from the Sun to the center of the galaxy is 
3.9^12 X 10~^, which is based on the 13 lensed giants 
that are concentrated near the galactic center. An 
overall optical depth, based on the entire data set, 
was estimated to be 2.43^q'45 ^ 10~^, where the lat- 
ter number is thought to be smaller because some of 
the main sequence sources were closer to us than the 
center of the galaxy. Both optical depths account for 
empirical determinations of the detection efficiency, 
which never exceeded 0.4 and decreased toward zero 
as event time scales approached a few days. These 
optical depths imply that far fewer than 100 million 
stars would have to be monitored to detect 350 events 
per year if the detection efficiency were increased by 
using a larger telescope at a better site along with 
only modest improvements in technology and obser- 
vational procedures (C. Stubbs, private communica- 
tion, 1996). If each event lasts an average of about 40 
days (Alcock, et al, 1996), 350 events in a 190 day ob- 
serving period would mean that almost 75 stars would 
be lensed at any one time. With continued improve- 
ments in technology and observing procedures, the 
observation of about 3000 microlensing events over 
about an 8 year period is not an unreasonable goal. 
These events would be covered with precise, high time 
resolution photometry about 58% of the time on av- 
erage. We will adopt these numbers in the examples 
to follow, where the number of detected small planets 
scales directly both with the total number of events 
and with the percentage of photometric coverage and 
the number of large planets with just the total num- 
ber of events. 

5. Sample expectations from a microlensing 
secirch 

What is expected from any microlensing search for 
planets depends entirely on what is assumed about 
the nature of likely planetary systems. It is guar- 
anteed that the distribution of planet properties we 
find is going to be unexpected, but to illustrate some 
examples, and to show the assumption dependence 
of the expectations, we shall base our model systems 
on our own solar system. In about half of the 3000 
assumed events, the lens will be a member of a bi- 
nary system although this binary nature will usually 



not be revealed. (Close binaries will act as a single 
lens and binaries separated by more than twice the 
Einstein ring radius of the more massive member will 
[usually] act as two separate lenses. Attempts to ob- 
tain a secure distribution of binary separations are 
frustrated by severe observational bias (D. Popper, 
private communication, 1996).) To be conservative, 
we assume that none of the binary systems have any 
planets, but to be equally optimistic we assume that 
all of the remaining 1500 events involve lenses with 
planetary systems similar to our own, but with only 
four or five planets. 

All 1500 remaining lenses have a Venus (to/Mq = 
2.5 X 10-6) at 0.7 AU, and an Earth (m/M© = 
3 X 10-6) at 1 AU, but only half or 750 lenses have 
a Jupiter (m/M© = 1 x lO^"*) at 5 AU and a Saturn 
(m/Mg = 3 X 10-'') at 10 AU. The other 750 lenses 
have Uranuses {tti/Mq = 5 x IQ-^) at the distances 
of the Jupiters and Saturns and, for these latter sys- 
tems, an extra Earth at 2.5 AU. It is difficult to make 
a Jupiter in the current planet formation paradigm 
within the observationally estimated lifetimes of pre- 
planetary disks, so we remove Jupiter from half of 
the systems. The extra Earth at 2.5 AU is motivated 
by Wetherill's (1994) calculations showing the per- 
sistence of Earth mass planets in the asteroid belt if 
Jupiter is absent. Given these systems around 1500 of 
the lenses, how many of the planets would be detected 
in the dedicated microlensing search described above? 
What is the nature of the minimum data set that 
would be obtained, and how well could it be inter- 
preted? Before attempting to answer these questions, 
we must make assumptions concerning the variation 
of the system properties with the mass of the central 
star. 

We shall determine the expected number of de- 
tected planets for three distinct sets of assumptions: 
1. The planet-star mass ratios and separations re- 
main invariant for all stellar masses M. This means 
that smaller stars have smaller planets, but that the 
Earths, for example, will always be at 1 (or 2.5) AU 
regardless of the value of M. There is some weak 
justification for assuming that the planet-star sepa- 
rations are the same regardless of the mass of the 
primary, at least for Jupiters, since the ice conden- 
sation point in a preplanetary nebula where Jupiters 
are thought to form may always be near 5 AU (Boss, 
1996). 2. The planet-star mass ratios remain invari- 
ant, but the separations are always the same fraction 
of the Einstein ring radius at z = 0.5. Here smaller 
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Fig. 7. — Distribution over the to, a plane of the prob- 
ability of detecting a planet P{m, a, Zg) during an on- 
going microlensing event for the case where the aver- 
age fractional distance to the lens Zs = 0.8 and where 
the planetary mass m and semimajor axis a are as- 
sumed invariant over the lens mass distribution. 

stars have smaller planets, but the planets move closer 
to their central stars as AI decreases. 3. The masses 
and separations remain invariant instead of the mass 
ratios and separations. In this last case, the Jupiters 
for example will always have m = O.OOIM0 for all val- 
ues of M and will always be 5 AU from their central 
stars. 

The results of applying Eq. for the three sets of 
assumptions and two values of Zg are given in Table 

and Fig. |^ shows the probability distribution over 
the TO, a plane for the model corresponding to the last 
column. In this application of Eq. (^, each planet 
type is treated independently, so the probabilities in- 
clude the perhaps unlikely situation of detecting two 
planets around the same star. Multiplying the num- 
ber of assumed planets by these probabilities yields 
the number of each of the planets that are detectable. 
We shall assume that the fraction of the detectable 
small planets (Earths and Venuses) that are actually 
characterized is the same as the 58% averaged maxi- 
mum coverage, but that all of the detectable Jupiters 
are detected with the probability curve in Fig. ^ now 
representing the average fraction of the more than one 



day long light curve perturbation that is followed with 
the high time resolution photometry. Saturn mass 
planets will also have perturbation time scales com- 
parable to a day and all will be assumed detected, but 
Uranus perturbations will have time scales compara- 
ble to the time the GC can be observed on an aver- 
age night at a single southern observatory, so perhaps 
75% of the detectable Uranuses may be found. These 
assumptions are reflected in the number of detected 
planets listed in Table ^. 

In the second of the three models in Table |^, where 
both the m/M and a/ Re are invariant for each of the 
planets as the lens mass varies, the probabilities sim- 
ply reflect the scaled values derived from Fig. H For 
example, the Jupiters in this model always correspond 
to a planet/lens mass ratio of 0.001 with a semimajor 
axis that is always about 1.3 times the Einstein ring 
radius when z — 0.5, where its probability of detec- 
tion is 0.175 for a lens of any mass. Hence, the integral 
over the mass distribution (Eq. (||)) yields just this 
probability in Table ||. The detection probability is 
reduced somewhat for z = 0.8 in the table because 
the semimajor axis is the same as for z — 0.5 but the 
lensing zone has moved in with the Einstein ring ra- 
dius leaving the Jupiters to the right of the peak in 
Fig. ^. Uranus at 5 AU for a solar mass star occupies 
the same position as Jupiter, but its probability of 
detection is scaled downwards by yJrajM to the ap- 
propriate value of 0.039 for z = 0.5 in the table. The 
probabilities for detecting Venuses and Earths nomi- 
nally at 0.7 and 1 AU respectively (Nominally because 
these are the semimajor axes only for A/ = 1 in this 
model.) are negligibly small because these planets are 
always too far inside the lensing zone in this model, 
although one of the Earths is detected at 1 AU for 
Zg = 0.8 and few more for both Zg = 0.5 and 0.8 at 
2.5 AU are detected. 

The semimajor axes are fixed in the other two mod- 
els, so the lensing zone scans the planetary positions 
as the lens mass is reduced. For the lower mass lenses 
this means that the Venuses and Earths at 0.7 and 
1 AU move into the lensing zone so more are de- 
tected. More are detected when the masses are held 
constant because the planet /lens mass ratio increases 
with decreasing lens mass. All of these very close 
planets would be detected about stars less massive 
than the Sun. The detection probabilities are higher 
for Zg = 0.8 compared to Zg = 0.5 for the planets that 
are close to their central stars and lower for those 
further away because moving the lensing zone closer 
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m,i/M^ a 


i ~const. 




TTii/M^ ail Re '^const. 


mi, 




-^const. 






Zs = 0.5 




Zs = 0.8 




Zs = 0.5 


Zs = 0.8 




Zs = 0.5 


Zs = 0.8 


planet 


Probability of detection 


Venus 


6.64 X 10" 


4 


1.45 X 10" 


'i 


6.84 X 10-^ 


1.30 X 10- 


4 


1.75 X 10" 


-3 


3.80 X 10- 


-3 


Earth(l) 


2.31 X 10" 


3 


3.75 X 10- 


-3 


2.08 X 10"* 


3.95 X 10- 


4 


5.92 X 10" 


-3 


8.98 X 10- 


-3 


Earth(2.5) 


7.65 X 10- 


3 


7.30 X 10- 


-3 


4.86 X 10-3 


8.37 X 10- 


3 


2.36 X 10- 


-2 


1.19 X 10- 


-2 


Uranus (5) 


1.78 X 10- 


2 


1.14 X 10- 


-2 


3.92 X 10-2 


3.00 X 10- 


2 


2.54 X 10- 


-2 


1.55 X 10- 


-2 


Uranus(lO) 


3.49 X 10- 


3 


1.68 X 10- 


-3 


9.80 X 10-3 


6.29 X 10- 


3 


4.30 X 10- 


-3 


1.89 X 10- 


-3 


Jupiter 


7.97 X 10- 


2 


5.10 X 10- 


-2 


1.75 X 10-1 


1.34 X 10- 


1 


1.14 X 10- 


-1 


6.95 X 10- 


-2 


Saturn 


8.65 X 10- 


3 


4.11 X 10- 


'3 


2.40 X 10-2 


1.54 X 10- 


2 


1.05 X 10- 


-2 


4.62 X 10- 


-3 



Table 2: Detection probabilities P(to, a, Zg) 





rrii/M, a 


i ^const. 


rrii/M, ai/Rs 


'^const. 




--•const. 


planet 


Zs = 0.5 


Zs = 0.8 


Zs = 0.5 




= 0.8 


Zs = 0.5 


Zs = 0.8 


Venus 


1 


2 










2 


3 


Earth(l) 


2 


3 







1 


5 


8 


Earth(2.5) 


3 


3 


2 




4 


6 


5 


Uranus(5) 


10 


6 


22 




17 


14 


9 


Uranus(lO) 


2 


1 


5 




4 


2 


1 


Jupiter 


60 


38 


131 




100 


85 


52 


Saturn 


6 


3 


18 




12 


8 


3 



Table 3: Number of planets detected. 



to the lens for the higher value of Zs moves the inner 
planets up the rising slope of the probability curve in- 
side the lensing zone, while moving the outer planets 
further down the decreasing curve outside the lensing 
zone. 

The entries into Tables || and ^ identify the plan- 
ets detected, but the data set obtained from the mi- 
crolensing survey yields only the planet-star mass ra- 
tio and the projected position of the planet in units 
of the Einstein ring radius from the light curve alone. 
How much information about the assumed distribu- 
tion of planetary masses and separations can be re- 
trieved from the actual data set? We begin to answer 
this question by constructing that data set from the 
distributions of the stellar masses and of the projected 
positions of the planets. Other information about the 
system that can be derived from different types of 
observations will be neglected for the time being. 

We shall develop the sample data set only for the 
model corresponding to the last column in Tables § 



and ^, where Zg = 0.8 and rrii and at do not change 
with M. To determine the likely planet/star mass ra- 
tios for the detected events, we note from Eq. that 
the fraction of detected planets of given assumed mass 
and semimajor axis that orbit stars between masses 
M = 0.08 and M is given by 

^ Cs^{M')Pim.,M',a,,Zs)dM' ^^^^ 
lo.os ^{M')P{m,,M', a„ Zs)dM' ' 

These cumulative fractions are shown in Fig. ^ as 
a function of M. These fractions depend only on 
the planet-lens separation and not on the planetary 
masses, so the curves are labeled only with the semi- 
major axes. As mentioned earlier, the Venuses at 
0.7 AU and the Earths at 1 AU are mostly detected 
around spectral class M stars, but the planets at 2.5 
AU and beyond are distributed over a wider range 
of lens masses. To assign the lens mass for each of 
the detected planets, we spread the detected masses 
over the distributions shown in Fig |^ in proportion 
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Table 4: Probable distributions of planet/star mass ratios and separations 
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Fig. 8. — Fraction of planets detected in the mass 
range 0.08 to M for each semimajor axis in the model 
planetary systems. Close planets are preferentially 
detected around the lower mass stars, whereas, the 
more distant planets are detected over a wider range 
of masses. 

to the fraction detected in each range. For exam- 
ple, from Table |3| the last column shows that about 
three Venuses would have been detected. From Fig. 
^ about 1/3 of the detected planets should be found 
about lenses in the mass range 0.08 < M < 0.10, 
another third in 0.10 < ill < 0.14 and the remain- 
ing third in M > 0.14. We can assign the probable 
lens masses for the Venus detections such that there 
is one within each range of masses corresponding to 
1/3 of the probable detections. If we choose these 
to correspond to values of the cumulative fraction of 
detections near 1/6, 1/2, 5/6, the lens masses about 
which the Venuses are detected would be 0.09, 0.12, 
0.20 respectively from Fig. ||. We can systematically 
assign lens masses to N detected planets with a par- 
ticular semimajor axis by choosing those values of M 
where the cumulative fraction of probable detections 
in Fig. I is 1/{2N) + k/N, fc = 0, 1, • • • iV - 1 as was 
done for Venus. The results of this exercise are shown 
in Table |. 

The determination of the distribution of normal- 
ized projected separations is more involved and not 
so well constrained as that for the planet/star mass 



ratios. An Earth with a = 2.5AC/ at its extreme pro- 
jected separation would be inside for a large star, 
but outside for a small star. But this Earth, uni- 
formly distributed on a sphere of radius o, could be 
at any projected separation from the lens less than a 
with a probability distribution given by p{r)dr in Eq. 

. For a particular lens mass M with a planet at a 
particular semimajor axis a, the probability of detect- 
ing the planet in a range Ar about r < a is the prob- 
ability of finding the planet in this range weighted by 
the probability of detecting the planet at that partic- 
ular separation. The first probability peaks at r = a 
and the second peaks near the Einstein ring radius for 
the point sources considered here. Although the prob- 
ability distribution P(m, a, M, Zs) obtained by scal- 
ing and shifting the curve in Fig. ^ has been aver- 
aged over the line of sight for a model distribution of 
stars, we can use it as an approximation to the sec- 
ond probability by replacing a hy r < a. It peaks 
near the Einstein ring radius and should give a rea- 
sonably good representation of the proper weighting 
of the projected distance probability. The fraction fd 
of planets of mass m, semimajor axis a that would be 
detected in the range 0— r about an ensemble of lenses 
of mass M, average fractional distance Zs is thus 

j^f^^^ ^ lo Pir')Pim,r',M, z,)dr' ^^^^ 
p{r')P{m, r',M, Zs)dr' 

In Fig. 1^ fdir) is shown for intermediate values 
of the lens mass M in the lens mass distribution as- 
signed to each planet type. Three Venus mass plan- 
ets are detected in our model data set about lenses 
with masses 0.09, 0.12 and O.2OM0 with correspond- 
ing Einstein ring radii of 0.97, 1.12, 1.44 AU respec- 
tively calculated at the average distance z^ — 0.8. 
The curves for the extreme values of M for Venus 
in Fig. P almost coincide with that for the interme- 
diate value shown for M = 0.12. We therefore dis- 
tribute the distances of the detected Venuses accord- 
ing to this curve (like the distributions of lens masses) 
at values of r/a = 0.82, 0.95, 0.98 corresponding to 
fd{r) ~ 0.16, 0.5, 0.83 respectively. We assign the 
closer planets to the smaller lenses since the closer 
proximity of the Einstein ring radii increases their 
probability of detection. The values of r/a and r/R^ 
are given in Table ^ 

The other planet types are treated in the same way 
with the separations being distributed according to 
fd{T) for a lens mass giving a curve approximately 
midway between the curves for the extreme values 
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Fig. 9. — Fraction of planets detected in range to 
r/a for each planetary type in the model systems. 
Curves for each planet type correspond to a lens mass 
near the mean value of the distribution of lens masses 
about which that planet type were found in Table ^. 
They are the basis for the assigned distribution in 
projected separations of planets. 

of the lens mass in the distribution. Here again the 
smallest separations are associated with the smallest 
lenses with the values given in Table |. As the semi- 
major axis gets small compared to the Einstein ring 
radius the distribution is dominated by the probabil- 
ity of the projected separation of the planet which 
peaks at r = a. The effect of the higher probability 
of detection at the Einstein ring radius becomes im- 
portant as the projected separation approaches this 
radius and beyond. Hence, we find that fd{r) for all 
the planet types moves toward the Venus curve when 
the lens mass is large. For the Earths at 1 AU at 
the largest lens mass M = 0.40, about which one of 
these closer Earths was detected, fdir) almost coin- 
cides with the Venus curve. For the Earths at 1 AU, 
the values of r/a and v/Re (at Zs = 0.8) in Table ^ 
are distributed according to the curve for M = 0.14. 
Similarly, the Earths at 2.5 AU are assigned sepa- 
rations according to the distribution for M — 0.25, 
Jupiters according to the distribution for M — 0.44, 
Saturns according to the distribution for M — 0.74, 
and Uranuses at 5 AU according to the distribution 



for M ~ 0.45. The single Uranus that was detected 
at 10 AU is placed at rja = 0.72 corresponding to 
/d(^) = 0.5 for M = 0.74. This procedure gives a 
reasonably probable distribution of separations, since 
the larger fraction detected at smaller distances for 
smaller lenses is accounted for by assigning the closer 
planets to the smaller lens masses. 

The values of m/M and r / Re in Table ^ comprise 
the information about the detected planets that can 
be learned from the light curve. Although additional 
information can be obtained by further observations 
of a different type, we start by seeing what can be 
deduced about the planetary systems from this infor- 
mation alone. 

6. Interpretation and discussion 

One of the first things that one notices about the 
data set in Table || is that the mass ratios for the de- 
tected planets are reasonably good signatures for the 
actual masses of the planets in spite of the varying 
lens masses. This results from the fact that the lens 
masses only span an order of magnitude whereas the 
planetary masses span three orders of magnitude. We 
would thus conclude that statistically m/M ^ 10~^ 
corresponds to terrestrial mass planets, mlM ~ 10^** 
to Uranus mass planets and m/M ^ 10^^ to Jupiter 
mass planets, where the terms terrestrial, Uranus and 
Jupiter are meant to constrain the masses only to 
about an order of magnitude from the least to the 
most massive in each designation. This grouping is 
evident in Fig. |l^ where we have plotted the val- 
ues of m/M vs t/Re- We also notice that the lower 
mass planets tend to fall inside the Einstein ring radii 
of their respective lenses, whereas the higher mass 
planets tend to fall outside. Thus, this solar system 
characteristic of our original assumptions is recov- 
ered. The concentration of the lower mass planets 
near t/Re ~ 1 are the Earths at 2.5 AU, and we 
could infer this approximate location for these plan- 
ets from the fact that the probability of detection is 
highest when the planets have semimajor axes near 
the Einstein ring radius and that radius is near 2.5 
AU for the dominant number of small mass stars in 
the distribution. The concentration of the medium 
and high mass planets near t/Re = 1.5 is also strik- 
ing, but this reflects the fact that these planets have 
semimajor axes that are larger than the Einstein rings 
of nearly all likely lenses, but have a greater proba- 
bility of being detected if their projected separations 
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Fig. 10. — Planet/lens mass ratios as a function of 
the projected separation for the constructed data set 
in Table |. 

are closer to the Einstein ring radius. Two Saturns 
and the single detected Uranus at 10 AU comprise the 
three points with t/Re > 2. 

The facts that all the sources will not be 8 kpc from 
the Sun and that the masses of real planets will not 
be limited to a few discrete values would introduce 
additional scatter into the data set displayed in Fig. 
p!o| . The groups would no longer be evident, although 
the general trend of the distribution of points would 
still prevail if other planetary systems were to have 
small planets close to their stars and larger planets 
further away. The very close large planets around 51 
Peg and 55 Cnc within .05 to .1 AU could not be 
detected, but the very large planet about 47 UMa at 
2.1 AU, if detected, would yield a point in Fig. |l^ 
most likely between values of t/Re between 0.4 and 
0.5 and m/M > 2.5 x IQ-^ (msini = 2.5x IO-^Mq). 
If many points were found in this region of Fig. |l^, 
it would imply the common existence of large planets 
considerably closer to their central stars than those 
in the solar system. 

The detection of sixteen planets with planet/lens 
mass ratios near 10~^ would imply that such plan- 
ets are relatively common given the low probability 
of their detection. However, care must be exercised 
in inferring the fraction of lenses having planets of 



various masses from the number of planets detected 
and the probabilities of detection, since the latter are 
dependent on the assumptions about the planetary 
distributions and the lens mass function. For exam- 
ple, with the same mass function, we detected respec- 
tively 8, 5 and 16 terrestrial type planets for Zg = 0.8 
for the three models, implying that our estimate of 
the fraction of lenses having terrestrial mass planets 
could not be constrained by less than a factor of 2. 
Even though a factor of 2 may not sound too bad, the 
small number of terrestrial planet detections makes 
such statistics suspect. In any case the model depen- 
dence of the interpretation is just as severe as that 
in deriving the expectations from the search. All the 
terrestrial planets are detected about stars consider- 
ably less massive than the Sun except higher mass 
lenses can be inferred for the terrestrial planets de- 
tected with r/i?E near 1. 

The determination of the fraction of lenses having 
Jupiters is similarly uncertain by at least a factor of 
2. However, if nearly all of these would correspond to 
values of r / R^ > 1 as in our one data set, we could 
infer the Jupiters to be relatively far away from their 
stars consistent with formation near the ice conden- 
sation point in the pre-planetary nebula. This consis- 
tency with the current paradigm of planet formation 
in the solar system and the inferred non dominance 
of giant planet migration (Ward and Hourigan, 1989; 
Lin, et al. 1996; Ward, 1996) in planetary system his- 
tories allows us the further inference of the existence 
of several terrestrial type planets inside the Jupiter 
orbit, even though the latter were not directly de- 
tected. The massive planets close to their primaries 
can be detected by radial velocity and astrometric 
techniques on the relatively short time scales of the 
orbit periods as already demonstrated. Statistics for 
massive planets in Jupiter-like orbits would be de- 
layed in such surveys, but would start accumulating 
immediately in an intensive microlensing search. An 
occasional detection of a low mass planet would sup- 
port the inference that where there are Jupiters far 
out, terrestrial planets are in close. 

The detection of the nine Uranus mass planets just 
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outside the Einstein Ring radii in Table || and Fig, 
would imply that the Jupiters are not there for these 
systems. A further inference could be that a relatively 
larger number terrestrial type planets may be inside 
5 AU for those particular stars than exist in our own 
solar system (Wetherill, 1994). The ratio of the num- 
ber of Uranuses to the number of Jupiters in such a 
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data set would test the implication of current theories 
that Jupiters at distances where they could form in 
the nebula may be even more rare than we assumed 
at the outset. Radial velocity and astrometric sur- 
veys can determine in a relatively short time if many 
Jupiters have migrated to or are otherwise positioned 
close to their stars. 

We see that several conclusions can be drawn from 
the statistics of the planet/star mass ratios and pro- 
jected separations obtainable from the microlensing 
light curves alone. The fraction of the stars that 
have planets of various masses in their lensing zones 
can be estimated. However, the model dependence 
of the interpretation of these data limits the accu- 
racy of the estimates to no better than a factor of 
2. The planet/star mass ratios inferred from the data 
are reasonably good indicators of the actual planetary 
masses because of the limited range of likely stellar 
masses. A rough dependence of the planet-star sep- 
aration on the planetary mass may be inferred from 
a plot of r/i?E vs m/M (Fig. [l^ ). The number of 
Uranuses detected where Jupiters might be expected 
constrains the success rate for the Jupiter-forming 
process. 

Few of the terrestrial type planets that are detected 
will be close enough to their stars for habitability. 
Recent developments indicate that a way to infer the 
frequency of such planets may be the distribution of 
separations of the massive planets from their stars. 
About 1 in 30 stars in radial velicity surveys to date 
seem to have giant planetary mass objects close to 
their primaries (W. Cochran, private communication 
1996). There are several theoretical arguments show- 
ing how Jupiters could be removed from their region 
of formation and, in some cases, account for the close 
planets. These include torques from spiral density 
waves in a persistent nebula generated by fully formed 
planets (Lin et al, 1996), by planetary billiards among 
giant planets formed too close together for stabilty 
(Rasio and Ford, 1996, Weidenschilling and Marzari, 
1996) or by rapid migration of the solid body accre- 
tional cores toward the star before acquisition of the 
gaseous envelopes (Ward and Hourigan, 1989; Ward, 
1996). 

The observations and the theories that attempt to 
account for them imply that most of the material that 
would have gone into terrestrial type planets could 
have been eliminated through interactions with the 
giant planets or their cores that invaded the terrestrial 
planet territory. Therefore, the greatest contribution 



that a microlensing search for planets may make may 
be to determine where the giant or subgiant plan- 
ets are relative to their central stars. The absence 
of many Jupiters or Uranuses detected outside the 
Einstein ring radii of their central stars may indi- 
cate that terrestrial-type planets in habitable zones 
are infrequent. On the other hand, finding most of a 
reasonably large number of massive planets relatively 
far from their stars would indicate that the scattering 
or migration mechanisms did not dominate planetary 
system histories and may give a better inference of 
the frequency of occurence of terrestrial planets than 
that provided by the meagre statistics of the latter's 
direct detection. 

We have limited the discussion so far to just what 
may be learned from the statistics of the microlensing 
light curves, since that is the only information that 
would be systematically obtained in the microlens- 
ing search described above. More information can be 
obtained about an individual system with additional 
observations. If the lens can be spectrally classified, 
the actual masses of the planets of such lenses can 
be obtained. Although spectral classification allows 
an estimate of the distance Dq^ to the lens from the 
known luminosity and estimates of interstellar extinc- 
tion, the additional unknown Dqs in limits the ac- 
curacy with which the actual separation between lens 
and planet can be determined. An isolated star with 
visual magnitude m„ = 24 can be spectrally classified 
with a Keck telescope (R.M. Rich, private communi- 
cation, 1996). In the crowded field toward the galactic 
center we should at least be able to classify a star with 
my = 22, which corresponds to a K5V main sequence 
dwarf at the galactic center. The time after the event 
that one must wait in order to separate the light of 
the lens from that of the source depends on the rate 
of development of interferometric techniques, but we 
would know the masses of those planetary compan- 
ions of lenses earlier than K5V. The mass of a K5V 
star is about O.69M0 (Allen, 1973), so the remaining 
lenses of later spectral type have a masses ranging 
over a factor of 8. 

This is not that great an advantage since only 
about 25% of the lenses in Table | have M > 0.69. 
Still, masses of nearly all of the detected planets 
would be known to within about a factor of 3 (an or- 
der of magnitude from the least to the highest mass 
estimate) from a microlensing search, and the pro- 
jected separations in terms of the Einstein ring radii 
would yield estimates of the actual separations. In 
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some circumstances observations using effects of dif- 
ferential limb darkening of finite sized sources (Gould 
and Welch, 1996) or of limb brightened H-K lines of 
Calcium (Loeb and Sassalov, 1995) as a small plane- 
tary Einstein ring traversed the source might be used 
to gain additional information about the lens system. 
However, discussion of the practical aspects of in- 
corporating such observations into an intensive mi- 
crolensing survey should await inclusion of the effects 
of finite sized sources in the detection of small plan- 
ets in an analysis like that above. Inclusion of the 
finite source size effects in a development similar to 
the above must await the accumulation of an array of 
detection probability curves like those of Bennett and 
Rhie (1996) covering a range of parameters. A pos- 
sible approach to small planet detection with finite 
sized sources is outlined in the next section. 

7. Finite sized sources 

Allowing for the finite size of the source changes 
the procedure used to determine the planet/lens mass 
ratio for small planets, since the time scale of the 
planetary perturbation is now extended beyond the 
transit time of a point source across the Einstein ring 
of the planet by the finite angular size of the source 
star. The maximum amplitude of the perturbation of 
the light curve is also reduced since a small planet can 
amplify only a fraction of an extended image. The de- 
termination of the probability of a planetary detection 
(Bennett and Rhie, 1996) during an event is much dif- 
ferent from that developed analytically by Gould and 
Loeb (1992). The Bennett and Rhie technique will be 
described below, and an example calculation of aver- 
aging a Bennett and Rhie detection probability over 
the stellar mass function establishes a procedure for 
determining the overall probability of detecting small 
mass planets. Even though this more correctly de- 
termined detection probability for small mass planets 
has little to do with that determined by scaling the 
Gould and Loeb probability curve for point sources, 
their remarkably close agreement for this one exam- 
ple implies that our conclusions above for the overall 
probability of detecting small mass planets may not 
be in error by more than a factor of 2. 

As pointed out above, the finite size of the source 
star must be taken into account when the Einstein 
ring angular radius of the planet becomes compara- 
ble with the angular size of the source star. The an- 
gular size of a star at 8 kpc is 0* = 6 x lO~^i?*/i?0 
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Fig. 11. — Probability of detecting a planet as a 
function of the projected lens-planet separation when 
to/A/ = 10^^ and the source star at the center of 
the galaxy has a radius of 3Rq. (after Bennett and 
Rhie, 1996). The Bennett and Rhie curves have been 
extended beyond r — 2Re by determining the proba- 
bility of detecting a distant, essentially isolated planet 
given that the planet's central star has acted as a lens. 



arcsec, whereas the Einstein ring radius of a planet 
for a source at 8 kpc is 1 x 10~'^-y/m(l — z)/{Mqz) = 
2 X 10~^ arcsec for a Jupiter at 6 kpc and 1 x 10^® arc- 
sec for an Earth at 6 kpc. So only for the largest super 
giants will the perturbation by a Jupiter mass planet 
show the effects of the finite size of the source, but es- 
sentially every perturbation by an Earth mass planet 
will be greatly affected by the source size. Bennett 
and Rhie (1996) determine detection probabilities for 
planet-star mass ratios of 10~^ and 10"'* for values 
of z of 0.5 and 0.8 with source radii of i?* = 3 and 
13i?0. We show how these probabilities can be aver- 
aged over the mass function, but that the probability 
determinations must be extended over the parameter 
space before rigorous overall probabilities of detection 
of Earth mass planets can be found for the finite sized 
sources. 

Fig. |l| shows the Bennett and Rhie (1996) prob- 
abilities of detecting a planet with m/M — 10^^ (an 
Earth mass when M = 0.3) for the two lens distances 
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as a function of the planet's projected separation from 
its central star when the source star at the center 
of the galaxy has a radius of 3Rq. This radius is 
a little larger than the peak in the distribution of 
radii of main sequence turn-off stars in the galactic 
bulge (Loeb and Sasselov, 1996). The probabilities 
are determined by numerically calculating many light 
curves for a set of impact parameters of the source 
relative to the lens spanning the Einstein ring of the 
lens and for all orientations of the source trajectory 
across the ring relative to the lens-planet line. The 
fraction of the light curves that displayed a perturba- 
tion more than 4% of the unperturbed amplification 
for more than 1/400 of the time scale of the entire 
event gives the probability of detection. The effect 
of the finite sized source is shown dramatically by 
the strong dip in the detection probability when the 
projected separation of the planet is near the Ein- 
stein ring radius. A point source yields a maximum 
detection probability of detection when the planet is 
very near the Einstein ring radius. This behavior for 
the finite sized source follows because the image is 
spread out to its maximum size when it is near the 
Einstein ring radius and the planet can focus only 
a much smaller fraction of the image light — failing to 
reach threshold detection. It is of interest to note that 
the maximum probability of detection for z = 0.5 is 
approximately double the maximum scaled probabil- 
ity for the point source case of Gould and Loeb (1992) 
(0.034 vs. 0.017 for m/M = 10"^), although the lat- 
ter would be increased somewhat if 4% instead of 5% 
were used as the detection criterion. The Bennett and 
Rhie detection probability falls somewhat for z = 0.8 
because the Einstein ring radius of the planet is re- 
duced at the larger value of z so a smaller fraction of 
the source image is amplified by the planet for similar 
geometries. 

Below we shall average these probabilities over the 
mass function of the lenses for planet orbital semi- 
major axes 0.7, 1.0, 2.5 and 5 AU respectively. The 
largest semimajor axis corresponds to a separation of 
5.5i?E for M = 0.08 and z = 0.8. We have therefore 
extended the Bennett and Rhie curves from r = 2Re 
to 5.5Re by estimating the probabilities of detecting 
the distant planet (Appendix B). 

The functional dependence of the probability curves 
is P = P(to, M, i?*, z, r), where r is the projected sep- 
aration of the planet from the lens in the lens plane 
and i?* is the radius of the source (assumed to be at 8 
kpc distance). The probability of detecting a planet 



of mass m with semimajor axis a during an event with 
source radius i?* , lens of mass M located at fractional 
distance z is the probability of finding the planet at 
projected distance between r and r + dr (Eq. (||)) 
times P{m, M, i?,, z, r) integrated from to a or 



P(m, M, R„,z,a) = / 
Jo 



" P{m,M, R^,z,r)rdr 



(12) 



1 - 



where all distances are determined in terms of the 
Einstein ring radius. So in general the integral is a 
function of M for given m, P*, z, a. If the probability 
curves are left in the broken line form shown in Fig. 
|lT| , each linear segment substituted into Eq. ( [l^ ) is 
analytically integrable and the singularity in the in- 
tegrand at r = a is removed. Similar to the Gould 
and Loeb probabilities, the probability of detecting a 
planet of mass m with semimajor axis a as a compan- 
ion of a lens at z in an event with source radius P, 
averaged over the event mass distribution is 

P(m, P*, z, a) = [ T{M)P{m, M, R^,z, a)dM, 
Jo. 08 

(13) 

The curves in Fig. |l| correspond to a particular 
mass ratio, so we can only evaluate the averaged prob- 
ability in Eq. ( p^ for the case where the planet/star 
mass ratio is fixed at this particular value as the lens 
mass varies. These probabilities are shown in Ta- 
ble 1^ along with those for point sources. With the 
exception of a = Q.7AU, the probabilities for point 
and finite sources are remarkably close. However, we 
cannot draw any inferences from these comparisons 
except perhaps to think that the probabilities for de- 
tecting Earth mass planets when point sources are 
assumed are closer to the real values than we might 
have expected. The finite source probabilities must 
be averaged over the distribution of lenses along the 
line of sight and over the distribution of sizes of those 
sources that will be monitored in any microlensing 
survey. This requires the generation of curves like 
those in Fig. |ll| for a sufficient number of source 
sizes and lens distances for reasonable interpolation 
over the distributions. In addition, if we were to con- 
sider the model where the mass of the planet is held 
constant instead of the mass ratio, detection prob- 
ability curves must be generated for various mass 
ratios. Hence, more accurate estimates of the de- 
tection of Earth mass planets through microlensing 
must await generation of this array of Bennett and 
Rhie type probability curves or, as may be unlikely. 
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P(TO,i?*,z,a) (Eq. (13)) 


P{m,a,Za)A (Eq. (9)) 


a AU 


z=0.5 


z=0.8 


z=0.5 


z=0.8 


0.7 
1.0 
2.5 
5.0 


0.000155 
0.00239 
0.0158 
0.0105 


0.000992 
0.00405 
0.0121 
0.00592 


0.00133 
0.00421 
0.0140 
0.00797 


0.00289 
0.00685 
0.0133 
0.00510 



Table 5: Probabilities of detecting a planet with m/M =10 ^ and i?, = 3i?0 with probabilities derived for point 
sources with m/M fixed at 10~^ included for comparison 



a clever scheme of avoiding the lengthy calculations 
can be contrived. On the other hand, the comparison 
of the probabilities in Table H with those from the 
point source calculations, may indicate that our esti- 
mates in Tables |^ and ^ for the Earth mass planets 
are correct to within a factor of 2. 

For the point source, the planet only perturbs the 
single lens light curve significantly when it is within 
a planetary Einstein ring radius of one of the unper- 
turbed image positions in the lens plane. If the angu- 
lar size of this Einstein ring is small compared with 
that of the source star, it will be small compared with 
the size of the image in the lens plane. Nothing much 
happens as the planetary Einstein ring approaches the 
edge of a large image, since it can amplify only a small 
fraction of the image light. That fraction of course 
depends on the angular size of the image (or angu- 
lar size of the source star) compared to that of the 
planetary Einstein ring, so bigger perturbations occur 
for smaller sources. Generally, a significant perturba- 
tion occurs only as the planetary Einstein ring moves 
well into the finite image, which occurs, at least for 
a planet outside the lens Einstein ring radius, as the 
caustic curve in the source plane moves across the face 
of the source star. The picture is more complicated 
if the planet is inside the lens Einstein ring radius 
since there are two caustics giving positive amplifica- 
tion separated by a deep negative trough, If the finite 
sized source covers both caustic and negative trough, 
the resulting cancellation can keep the total amplifi- 
cation undetectably small (Bennett and Rhie, 1996). 
If the source angular diameter is less than the sepa- 
ration of the two caustics a negative perturbation of 
the light curve can be detected for the inside planet. 

In most cases the perturbation of the light curve 
will have more structure than that shown by the ex- 
ample of Fig. H, since regions of positive amplifica- 



tion are always interleaved with regions of negative 
amplification. Examples of such light curves for sev- 
eral planet/lens mass ratios, several geometries and 
several source sizes are shown by Bennett and Rhie 
(1996) and by Wambsganss (1996). The perturba- 
tion light curve structure is most pronounced when 
the source size is smaller than the caustic in the 
source plane, but persists if the source is not too 
large, where a typical pattern is a negative perturba- 
tion surrounded by smaller positive perturbations for 
a planet inside the Einstein ring of the lensing star and 
a positive perturbation surrounded by smaller nega- 
tive perturbations for a planet outside the Einstein 
ring (Bennett and Rhie, 1996). 

The structure of the light curve perturbation com- 
plicates the simple picture of determining the planet/ 
lens mass ratio from the ratio of the perturbation time 
scale and the event time scale, but at the same time 
supplies sufficient constraints for a more precise deter- 
mination of that mass ratio. The normalized impact 
parameter of the lens-source encounter is known from 
the maximum of the light curve, and the angle at 
which the source encounters the lens-planet line and 
the location of the unperturbed images of the source 
are known from the position of the perturbation on 
the light curve. The projected position of the planet 
must be near one of the unperturbed images in order 
effect the perturbation. The selection of the correct 
image is usually very easy with a well sampled light 
curve as noted above (Gaudi and Gould, 1996). The 
planet may pass either side of the image and create 
the same maximum perturbation, but with different 
symmetries in the wings of the light curve that are 
nearly always discernable for passage near the image 
inside the Einstein ring (minor image) and usually 
discernable for passage near the outside image (ma- 
jor image) with 1% photometry and good coverage 
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(Gaudi and Gould, 1996). If the ratio of the source 
angular radius to the angular radius of the planetary 
Einstein ring is not too large, a condition satisfied 
by Earth mass planets when the source is a turn- 
off main sequence star at the center of the galaxy 
(i?* ~ 2to3i?0), the remaining uncertainty in the 
projected position of the planet is negligibly small. 
As the angular size of the source is known spectro- 
scopically from the apparent brightness and an esti- 
mate of interstellar absorption, the planet/star mass 
ratio, as the remaining free parameter, can be ad- 
justed to match the light curve and thereby constrain 
its value quite well. A model of limb darkening and 
slight adjustments in the source trajectory through 
the lens-planet system may be necessary to refine the 
determination. 

About 25% of the events toward the galactic cen- 
ter will involve so-called clump giants (i?* ~ 13i?o) 
as sources (Alcock, et al 1996). Gaudi and Gould 
(1996) call attention to a potentially serious uncer- 
tainty in the planet/star mass ratio by demonstrat- 
ing a set of similar, low amplitude perturbation light 
curves involving the major image with the same maxi- 
mum change in amplification and the same full-width- 
half-maximum of the main part of the perturbation 
but corresponding to an order of magnitude range 
of planet /star mass ratios. The uncertainties are re- 
duced to acceptability with good coverage of both 
wings of the light curves with 1% photometry. How- 
ever, the practical difficulties in getting such cover- 
age routinely means that the uncertainties would re- 
main for many of these events. One means of reduc- 
ing the uncertainties in the mass ratio to acceptable 
values offered by Gaudi and Gould is to distinguish 
the curves with simultaneous photometry in the in- 
frared and optical. Giant stars are less limb darkened 
in the infrared than in the optical, and the distinct 
color effects of the event can select the proper mass 
ratio. On the other hand, the planetary perturba- 
tions are more difficult to detect in the first place 
with giant sources because of the lower amplitude of 
the perturbations. In most of the cases were the un- 
certainty would be severe, the observing strategy used 
by Bennett and Rhie (1996) would have failed to de- 
tect the planet (D. Bennett, private communication 
1996). This would eliminate most of the cases where 
the uncertainty in m/M was excessive from the data 
set. But even with those that remain, there are ways 
to determine both m/M and the projected separation 
to sufficient accuracy (Gaudi and Gould, 1996). 



So the effect of finite source size will not hinder the 
determination of the planet/lens mass ratio or the 
normalized projected separation of the planet from 
the lens. It thus makes sense to pursue the exercise 
outlined above to determine the overall probability of 
detection of small mass planets when the finite source 
size is accounted for. In the meantime, we can take 
some assurance from the comparison of the detection 
probabilities for finite and point sources in Table |, 
that our estimates for the detection of small planets 
may be within a factor of 2 of those that will be de- 
rived for the finite sources. 

We end by cautioning the reader to keep in mind 
the severe assumption dependence of the expectations 
for planet discoveries in any microlensing search. (Ta- 
ble H gives approximately 56, 138, or 81 planets de- 
tected for the three sets of assumptions assuming that 
no planets are found about members of binary star 
systems.) But he or she should also keep in mind that 
the detection scheme is robust with simple, ground 
based technology, and that if they are there, the plan- 
ets will be detected with useful bounds placed on their 
masses and planet-star separations with a data set 
based on the light curves alone. Finally, although 
microlensing is the only ground based scheme that 
is sensitive to Earth-mass planets, finding most of a 
large number of giant or semi-giant planets beyond 
the Einstein ring radii of their stars may be a better 
indicator of the frequency of occurrence of terrestrial 
planets than the few of the latter that are directly 
detected. 

Appendices 



A. Analytic Representation of the Gould and 
Loeb (1992) detection probability curve 



0.25a 0.25 



M V ^(1 - 



0, 0<j:<0.1 



= 0.001763066 



,7.6653342*(£c-0.1) 



- 1.0 



O.K X < 0.564 



= 0.4733727a; - 0.2069822, 0.564 < a; < 0.733 
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-84.9463918a;* + 286.1509291x3 - 361.2286639x2 
+202.6160685a; - 42.4667434, 0.733 < x < 0.9 



0.0475524d0 * x + .1204028, 0.9 < x < 1.043 

-1.3858063a;* + 6.0121987^3 - 9.8567924^^ 
+7.2772186a; - 1.8790527, 1.043 < a; < 1.4 

-0.1806818a; + 0.4164545, 1.4 < x < 1.664 

0.02869257X* - 0.3053152x3 + 1.23616560x2 
-2.2872802X + 1.6857629, 1.164 < x < 3.0 

-0.015d0*x + 0.075, 3.0 < x < 5.0 

0, X > 5.0 



B. Probability of detection of a planet that 
is far from its central star 

Bennett and Rhie (1996) require an amplification 
of 1.58 before a microlensing event is to be followed 
with high time resolution photometry in search of a 
planet. This requires that the source have an impact 
parameter relative to the lens of less than 0.7637 ac- 
cording to Eq. (^). A distant planet (a >> Re) will 
act as a solitary lens for most of its probable projected 
distances from the lens. The planetary Einstein ring 
radius Re is small compared to the projected ra- 
dius of the source star, which requires that the planet 
actually traverse the face of the source star for signif- 
icant amplification. (This compares with its having 
to traverse the image of the source star for signifi- 
cant amplification when it is close to its central star.) 
The maximum amplification is found by aligning the 
planet with the center of the source, multiplying Eq. 
(H) by an element of area of the source (normalized 
by the Einstein ring radius) integrating over the to- 
tal area of the source and dividing by the area. This 
yields 



A — 



402 

l<Ep 



where 9^^^^ is the angle subtended by the Einstein ring 
radius of the planet and 0^, is that subtended by the 
source radius i?*. For the mass 10~^ x O.O8M0 with 
i?* = 3i?0, Ajnax = 1-60 for z = 0.5, and A^ax = 
1.125 for z = 0.8. In both cases the planet would be 
detected by the 4% criterion, but only after the event 
because of the A = 1.58 threshold before a planet 
is sought. These maximum amplifications also show 
why it is necessary for the planet to be on the disk of 
the source star before detection is possible. 

The planet must be within the swath that is 2 x 
.7637i?i3 wide and 2a long centered on the lens star 
to have a chance of being intercepted by the source 
during the event. However, since the distant planet 
can be detected only after the event, we only consider 
the probability Pg that the planet is in one half of this 
swath. 



P. 



r dr 



a?^l — r"^ / a? 



dr 



dO 



^R/r a'^y/l 



P, = 0.0494 for z = 0.8 where the swath is 2 x 0.915 x 
5 (AU)2, and Ps = 0.0635 for z = 0.5 where the swath 
is 2 X 1.14 X 5(AU)2. 

For z — 0.8, 3i?o projected onto the lens plane 
and divided by Re for M = 0.08 is = 0.0122. 
If the planet is within the swath area it is approx- 
imately uniformly distributed in the direction per- 
pendicular to the long axis of the swath. The prob- 
ability that the source will encounter the planet is 
then r,/0.7637 = 0.016 and the probability of detec- 
tion is this probability times the probability that the 
planet is within the swath or 0.00079. This estab- 
lishes the point (5.46,0.00079) to which the Bennett 
and Rhie graph for z — 0.8 is extended in Fig. For 
z = 0.5, = 0.00614, the probability of interception 
is 0.00803 and probability of detection is 0.00051. So 
the Bennett and Rhie graph for z = 0.5 is extended 
to the point (4.37, 0.00051) in Fig. |ll|. 
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